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Abstract 

We construct fully backreacted holographic superfluid flow solutions in a five- dimensional 
theory that arises as a consistent truncation of low energy type IIB string theory. We 
construct a black hole with scalar and vector hair in this theory, and study the phase diagram. 
As expected, the superfluid phase ceases to exist for high enough superfluid velocity, but 
we show that the phase transition between normal and superfluid phases is always second 
order. We also analyze the zero temperature limit of these solutions. Interestingly, we find 
evidence that the emergent IR conformal symmetry of the zero-temperature domain wall is 
broken at high enough velocity. 



1 Introduction 

In the last few years there has been an intense effort to model superconductor/superffuid 
phase transitions using the AdS / CFT correspondence. The basic observation that makes this 
industry possible is the fact that at finite charge density and at sufficiently low temperatures, 
an AdS black hole in the presence of a charged scalar field is unstable to the formation of 
hair pp. Using the basic AdS/CFT dictionary [2j [3j S] , this gets easily interpreted as a 
superfluid-like phase transition in the dual field theory, cf. Weinberg [5]. 

Much of the work on holographic superconductors is done in the context of phenomeno- 
logical models, along the lines of the proposal originally presented in [6] . This is based on the 
minimal set-up of a charged massive scalar minimally coupled to Einstein- Maxwell theory. 
While many interesting results can be obtained within this minimal framework (see [H EJ [9] 
for reviews and references), such a bottom- up approach has some intrinsic limitations. Since 
the hope is that holographic constructions may eventually shed some light on some basic 
properties of high-T c superconductors, it would be desirable to have a microscopic under- 
standing of the underlying theory. This is something that phenomenological models, by 
definition, cannot offer. Secondly, they do not guarantee the existence of a quantum critical 
point in the phase diagram, which is instead expected to control the physics of high-T c su- 
perconductors. Indeed, the phenomenological models that one typically works with have no 
potentials but the mass term. However, it is expected that to have an emergent conformal 
symmetry in the infrared in the zero temperature limit, one should have potentials that allow 
symmetry-breaking minima [10]. Recently, some progress has been made in this respect and 
several microscopic embeddings of holographic superconductors have been proposed in the 
framework of type IIB string theory [TT], M-theory [T2], and D7-brane models [T5]. In these 
models, the potentials quite generically allow symmetry breaking vacua. 

Most studies have also been performed in the probe approximation, which is a large- 
charge limit in which the backreaction of the matter fields on the gravitational field is 
negligible. While many interesting results can be obtained with such a simplified setup when 
the temperatures are near the phase transition, the analysis becomes less and less reliable at 
very low temperatures, where the backreaction is non- negligible. This prevents exploration 
of interesting low temperature phenomena: in particular, understanding the ground state 
of holographic superconductors is outside the regime of applicability of the probe limit. 
Therefore, it is useful to realize holographic constructions where the backreaction is taken 
into account. Progress in this direction began with [14], where a (numerical) backreacted 
solution for the phenomenological model of [6] was presented. 

In trying to explore the phase diagram of holographic superconductors, an interesting 
direction was pursued in [T5J \W\ where the original holographic superfluid was studied in 
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the presence of a non- vanishing superfluid velocity (aka superfluid flow). Holographically 
this needs a non-trivial profile for a spatial component of the gauge field, besides the ever- 
present temporal component. The latter corresponds to a charge density and is necessary 
to have a phase transition in the first place (see [T7j for a recent alternative proposal). Two 
interesting results obtained in [15j [16] were to show the existence of a critical velocity above 
which the superfluid phase ceases to exist, as expected for physical superfluids, and the 
existence of a tricritical point in the velocity vs. temperature diagram where the order of 
the phase transition changes from second to first. Moreover, it was noticed in [18, 19] that 
these solutions can be efficiently compared to 2+1-dimensional superconducting thin films 
or wires. They behave very much like superfluids, in that an applied external magnetic 
field does not get expelled as if the gauge field were not dynamical. The four- dimensional 
gravitational model of [151 US] was further analyzed from this latter viewpoint in [19], where 
the system was in fact studied at fixed current rather than at fixed velocity. This choice 
allowed new checks, and remarkable agreement with some peculiar properties of real-life 
superconducting films (see |20j) was found. 

All solutions presented in [151 EH HH] have been obtained in the probe approximation. 
Hence, while being able to confront phenomena near or right below the critical temperature, 
not much could be said about the low temperature regime of such superfluid flows. This 
problem was addressed more recently in [21] , where the backreaction of the phenomenological 
four- dimensional model of [T5| [T6] was obtained. 

1.1 Summary of Results 

In this paper we take some concrete steps forward in the above program on superfluid flows: 
we focus our attention on models with known microscopic embedding and symmetry breaking 
vacua, and work at the backreacted level. Specifically, we will describe a holographic super- 
fluid flow in four dimensions by means of a fully backreacted solution of a five- dimensional 
gravitational system whose action arises as a consistent truncation of type IIB string theory 
[TT] . The effective theory is essentially Einstein-Maxwell theory with a Chern-Simons term, 
interacting with a complex charged scalar with a non-trivial potential. It can be obtained 
upon compactification of type IIB theory on an AdS^ x Y geometry, Y being a Sasaki-Einstein 
manifold. Using the numerical solutions that we find, we analyze several aspects of the rich 
phase diagram of this system. In particular, we present the plots of the scalar condensate 
against temperature and its dependence on the superfluid velocity, analyze the nature of the 
phase transition computing the free energy difference between the superconducting and the 
normal phase, and give some predictions on the zero temperature limit. 

As one would expect on physical grounds, we observe that for high enough velocity the 
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system stops superconducting. Interestingly, we find that for all velocities we have investi- 
gated the phase transition in these type IIB constructions is always second order. Hence, we 
do not find the tricritical point which characterizes the phase diagram of models with large 
charges. The same behavior was observed in the phenomenological but backreacted AdS^ 
model of [21] for low values of the scalar charge (in fact, only for the case q = 1 in their 
notation). The persistence of the second order phase transition has been observed also in 
the unbackreacted case for large masses of the scalar in five dimensions [22]. We will have 
some more comments on this in section HI 

One of the advantages of having a fully backreacted model is that one can also investigate 
the low temperature limit. In the zero velocity case, it is known [231 HH] that the type IIB 
hairy black hole solution tends to a domain wall with an emergent conformal symmetry in 
the deep IR. (This is in contrast with the phenomenological model of [24J where the potential 
has only a mass term and no symmetry breaking minima, and the zero temperature limit 
generically does not lead to an IR AdS geometry.) When the velocity is turned on and it 
is high enough, we find evidence that the solution stops being AdS in the IR. This suggests 
that beyond some critical velocity the IR conformality is lost. Along the way, we also discuss 
the importance of the frame comoving with the superfluid flow in these results. 

The rest of the paper is organized as follows. In section [2] we present the truncated type 
IIB five-dimensional action and the equations of motion. Our ansatz for the relevant fields, 
and the procedure we pursue to obtain numerical solutions with the desired features are 
discussed in section [3j Using these solutions, in section H] we study the phase diagram of the 
superfluid flow. In particular, we analyze the nature of the phase transition as a function of 
the superfluid velocity. In doing so, we compute the free energy of the superfluid phase and 
compare it to that of the normal phase (which is described, holographically, by a Reissner- 
Nordstrom black hole with no scalar hair). Finally, in section Owe study the T — > limit of 
some geometrical quantities like the Ricci scalar and the Riemann tensor squared. We also 
study the variation of the superfluid fraction as the temperature is lowered. These analyses 
allow us to explore the nature of the ground state of holographic type IIB superfluid flows. 
The appendices contain more technical material which might help the reader in following 
our analytical and numerical computations more closely. 

2 The IIB Set Up 

In [TJ] a consistent truncation of type IIB supergravity was presented, which has the 
structure of an Einstein-Maxwell (plus Chern- Simons) system in five dimensions coupled to 
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a charged scalar field with a non-trivial potential. The action reads 
Sun I d 5 Xy/^g 



I? 
3 



R - — F ab F ab + - I— j e abcde F ab F cd A e + 
-7i [ (dail>) 2 + sinh 2 - 2A a ) 2 - cosh 2 ( £ ) (5 - cosh^) 



• (2.1) 



Here, e 01234 = l/^/—g, and we have written the charged (complex) scalar by splitting the 
phase and the modulus in the form ipe l9 . For later convenience we recall that the Abelian 
gauge field A is dual to an i?-symmetry in the boundary field theory [11] and the scalar field 
has i?-charge R = 2. 

The matter equations of motion are 



^djh 2 ^F ab - ^L 3 ^e abcde F cd A e j 
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-.d a (V^d a iP) - - smh2i>(d b 9 - 2A b f + 
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^sinh-?/'(5 — cosh-?/') — 2 cosh 2 (^T^J sinh^ = 0. (2.3) 



The Einstein equations can be written as 

R-ab — ^9abR ~ ~L 2 (^F ac F b C — ^-F Cd F c dj + 



~2^ah+ ~ A 9ab^a ~ 7TF^9ab COsh 
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(5- cosh V>) =0,(2.4) 



with E ab = d a i)d b i) + smh 2 ^(d a 6 - 2A a )(d b 6 - 2A b ) . 



It is convenient to use the gauge invariance to shift away the angle 9 and also write the 
various expressions in terms of covariant derivatives. This basically means that we set 9 to 
zero in the above equations and use 

V a (^L 2 F\ - ^Lh\ cde F cd A e ) + ^L\ mr \F m F rs - 4 sinh 2 rj, A b = , (2.5) 

V a V>-2sinh2^(A 6 A fe ) + ^ (sinh ^(5 - cosh?/;) -2 cosh 2 \J?J sinh^) =0, (2.6) 

as the matter equations of motion. The leading terms in the scalar potential take the form 

12 3ip 2 
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+ ... 



(2.7) 



which have the immediate interpretation as the AdS cosmological constant and the scalar 
mass term. Typically, in a minimal phenomenological model the scalar potential has just 
the above two terms. Higher order terms affect mostly the very low temperature regime 
where the condensate becomes larger and thus the type IIB model can become substantially 
different from the minimal one. There are then two reasons as to why one should try and 
work out a fully backreacted solution for this type IIB model. The first is that the scalar has 
charge R = 2 and hence the probe approximation, which is a large charge scaling limit, is 
potentially inappropriate already at temperatures near the critical temperature. The second 
is that a backreacted solution would let one study the system in a regime (i.e. very low 
temperatures) where, as just noticed, the differences of the action (12.1 p with respect to that 
of a phenomenological model, are more apparent. 

Note that the scalar mass is m 2 = —3. In d — 4, this mass is in the range where 
the leading fall-off at the boundary, which is 0(l/r), corresponds to a non normalizable 
mode. So, using the AdS/CFT map, we will interpret it as the source of the dual field 
theory operator O. The subleading fall-off is 0(l/r 3 ) and corresponds to a condensate for 
O (whose dimension will therefore be A = 3). It is evident from the value of the R-charge 
and this fall-off that A = 3|i?|/2 and O is therefore a chiral primary [TTj . 

3 Hairy Black Hole Solution 

We want to construct a fully backreacted hairy black hole solution, holographically de- 
scribing a superfluid flow. To achieve this we must keep the metric (also) unfixed and find 
a self-consistent solution for the metric, the gauge field and the scalar. To have a charged 
scalar condense, we need to turn on both the scalar and the time component of the gauge 
field in the bulk [I]. Moreover, to obtain a non- vanishing superfluid flow, we should break 
the isotropy in the boundary directions that was present in the original holographic super- 
conductor construction of [TT] . Indeed, the superfluid velocity in (say) the x-direction is 
captured by the leading fall-off of the bulk gauge field component A x at the boundary, which 
should therefore have a non-trivial bulk profile. Altogether, this means that we need ip, A t 
and A x to be non-trivial. Since we would like to work with ordinary as opposed to partial dif- 
ferential equations, we look for an ansatz where these are functions purely of the holographic 
direction r : fortunately, this turns out to be enough to obtain a solution. Consistency of 
the Einstein equations then demands that we choose a metric ansatz of the form 

ds 2 = rt-t-dt 2 + i dr 2 - 2C{r)—dtdx + —B{r)dx 2 + — dy 2 + — dz 2 . (3.1) 

L z r z t(r) L z L z L z L z 
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The metric contains four independent functions, f(r),h(r),C(r) and B(r). Together with 
the ansatz for the gauge field and the scalar 

A = A t (r) dt + A x (r) dx , ip = ip(r) , (3.2) 

this will give rise to a set of seven independent equations for seven unknowns. Our ansatz 
here is essentially of the same form as the one in [21], albeit in one more dimension. This 
can be demonstrated by going over to an Eddington-Finkelstein form and working in a frame 
where the normal fluid considered in [21] is at rest. 

Let us first notice that with this choice of ansatz the terms in the equations of motion 
(12. 5p arising from the e abcde piece all vanish. A second important fact is that there are several 
scaling symmetries one should be aware of. In particular, the ambiguity in the units at the 
boundary for the time t and the distance along x translate to the following two scaling 
symmetries of the resulting equations 

t -> t/a , / -> a 2 / , /i4a/i, C^aC, A t ->■ ai f , (3.3) 

x x/b , B ->■ b 2 B , C -> b C , A, ->• b A x . (3.4) 

These are symmetries of the action and therefore of the equations of motion. Two further 
scaling symmetries of the system that we will use are 

(r,t,x,y,z,L) -> a(r,t,x,y,z,L) , {A t ,A x ) -> (A t ,A x )/a , (3.5) 
r->/3r, (t,x,y,z) ->• {t,x,y,z)/P , A^) /3(A t , A,) . (3.6) 

The first scaling changes the metric by a factor a 2 and leaves the gauge field invariant, but 
its effect is to scale the action (12. ip by an overall constant factor a 2 , therefore leaving the 
equations of motion unaffected. The second scaling is the usual holographic renormalization 
group operation in AdS, and it is easily seen that the metric, gauge field and the equations 
of motion are left invariant. Using the symmetries ( 13. 5 p and ( 13. 6 p we can scale the horizon 
radius r# and the AdS scale L to unity. We will assume this has been done in what follows, 
unless stated otherwise. 

The strategy we pursue to construct (numerically) our solution is as follows. First, using 
our ansatz, one can massage the equations of motion and end up with first order differential 
equations for / and h and second order differential equations for B, C, A t , A x and if). All in 
all we have then two first order and five second order equations resulting in twelve degrees 
of freedom. Therefore, to fix a solution we need twelve pieces of data. 

We start by considering the fields (I3.ip - (l3.2p near the horizon (r = r#) and expand their 
several components $ in a Taylor series as 

$ = $^ + $f(r- r// ) + ... . (3.7) 
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Requiring regularity of the solution at the horizon amounts to setting some specific coeffi- 
cients to zero. To linear order in (r — r#), the expansion at the horizon takes the form 



./ 


= f"{r-r H ) + ... 


(3.8) 


h 


= h% + hf(r-r H ) + ... 


(3.9) 


B 


= B* + B*(r-r H )+... 


(3.10) 


C 


= C*(r-r H ) + ... 


(3.11) 


A t 


= A* 1 {r-r H ) + ... 


(3.12) 




= A^ + A^(r-r H ) + ... 


(3.13) 


i> 


= ^ + <(r -r H ) + .... 


(3.14) 



That is, demanding regularity is tantamount to setting f^, C$ and (p^ to zero. Imposing 
now the equations of motion has the effect of putting further constraints on many coefficients, 
which all end up being determined by a small set of independent horizon data. It turns out 
that the coefficients can all be determined in terms of six independent data 

(h^B^C^A^A^). (3.15) 

This means that the solutions that we will find by integrating from the horizon will be a 
six-parameter family. All other coefficients are functions of these ones. One such relation 
which will be useful later is 

f? = (? + 2 C08h < - S^Ml) _ M . (3,6) 

The next step is to integrate the solution from the horizon out to the boundary (r — > oo), 
starting with the free horizon data f)3.15p . trying a suitable ansatz for the asymptotics of 
the fields at the boundary. In fact, the asymptotic expansion in five dimensions is subtle 
because, as already noticed, the mass of the scalar is such that there is a non-normalizable 
mode. To accommodate a generic solution obtained by integration from the horizon, we 
therefore need to turn on the non-normalizable mode of the scalar as well at the boundary. 
The non-normalizable mode triggers further logarithmic terms in the asymptotic expansion, 
so we need to keep track of them as well. It turns out that a combined series expansion in 
both l/r n and log r/r m 



$ = 

n=0 m=0 



works nicely. 
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Using a shooting technique we select, out of all possible solutions, those which match our 
physical requirements. In particular, we ask that the space be asymptotically AdS and that 
the source term for the field theory operator dual to the scalar field be vanishing, since we 
want the U(l) breaking to be spontaneous. 

We have found that the following asymptotic expansion solves the equations of motiorJl], 
while being general enough to match the curves arising from the integration from the horizon 

/ = tf + £ + £]°gr + ..., ^o + ^ + ^ + ^logr + ..., (3.18) 

B = B + ^ + ^lo g r + ..., C = C + ^ + ^lo g r + ..., (3.19) 

A t = A tfi + ^f + ^fhgr + ..., A x = A xfl + ^f + ^f\ogr + ..., (3.20) 

ip = h -7 + -^logr + ... . (3.21) 



Of course, not all of the above coefficients are independent. We relegate the explicit expres- 
sions for the dependent ones to Appendix |A] We merely note that when the non normalizable 
mode ipi is set to zero, the expressions are such that all the logarithmic pieces vanish as ex- 
pected. It can also be seen that the independent parameters at the boundary can be taken 
to be 

(ho, / 4 , £0, #4, C , C A , A t>0 , A t ,2, A xfi , A X)2 , ipi, . (3.22) 

To get asymptotically AdS solutions, we must set B , h to 1 and C , V>i to zero. The scaling 
symmetries can be used to accomplish the first two conditions, whereas we need to shoot for 
the last two. We are therefore left with eight independent boundary data. They are 

(/ 4 , B 4 , C 4 , A tfl , A t;2 , A x ,o, A x>2 , ip 3 ) . (3.23) 

We see here that the physical requirements we impose at the boundary do not fix as many 
integration constants of the ODE system, as the regularity conditions at the horizon does. 
This means that for the solutions that we obtain, there are hidden relations between the 
boundary data. Concretely, since there are only six independent pieces of horizon data this 
gives us relations between the above eight variables, which will then be used to study the 
phase diagram of the boundary theory. 

An important quantity for studying the thermodynamics of the system is of course the 
superfluid temperature. This corresponds to the black hole Hawking temperature, T. From 



1 What we do is to plug this expansion into the EoMs and demand that the result be zero order-by- 
order. We find that either this is satisfied identically or that the resulting relations can be interpreted as the 
definitions of higher order terms in the expansion. 
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the structure of the metric (13.11) we easily get 

T - r ^ ^ rg ) (3 24) 

1 AitI?h{r H )> {6 ' 

which is then also determined in terms of our horizon data. After some simple algebra, 
recalling we have set r# = 1 and using the horizon relation (13 . 16j) . we get 



T = — 

47T 



h 



H 



(3.25) 



4 Superfluid Flow Phase Transition 

We plot the result for the condensate versus the temperature in figure 1, for different 
values of the superfluid velocity. For the rest of the paper, we will introduce the notation 

where \x is the field theory chemical potential, O the (condensing) chiral primary operator, 
and £ the superfluid velocity in units of the chemical potential. When we work in an ensemble 
with fixed chemical potential, the meaningful (dimensionless) quantities relevant for the 
condensate plot are 

- and (4.2) 

In constructing the plots, we have also rescaled by the (velocity-dependent) factor yjl — £ 2 , 
which is nothing but the relativistic boost factor. 

From the form of the curves in figure 1, it is evident that there is a phase transition to 
a hairy black hole at low temperatures. As expected, the critical temperature decreases as 
the velocity is increased. For instance, for £ = 1/2 (which is the highest velocity we have 
investigated) we observe that T c (£ = 1/2) = 0.067 T c (£ = 0). It is clear from the condensate 
plot that the superfluid phase cannot exist for velocities that are much higher than this. 

One can compare the free energy of the normal phase (which corresponds to a Reissner- 
Nordstrom black hole with no hair) and the hairy/superfluid phase to see that the superfluid 
phase is favored when it exists. We collect some details of the free energy computation in 
appendix [Dj while figure 2 contains the free energy comparison between the superfluid phase 
and the normal phase at the same value of T In terms of S Iea defined in eq. ( IB. lip , the 
precise quantities we plot are 

<? O T 

vs. — . (4.3) 



^ 4 Vol 4 /i 4 ' ii 

The plot demonstrates that the phase transition stays second order for all values of the 
velocity, up to our numerical precision. This should be contrasted to the unbackreacted 
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Figure 1: Condensate plots for various values of the velocity £ = 0,0.1,0.33,0.4,0.5 (from 
right to left). The zero velocity case, £ = 0, which we report for ease of comparison, precisely 
agrees with existing results in the literature [TTj . 

cases previously considered in the literature, where the phase transition typically changes to 
first order for high enough values of the velocity [151 HSl [19] . In [2TJ a backreacted superfluid 
in AdSi was considered and it was found that for low enough values of the charge of the 
scalar field, the phase transition remained second order. Our type IIB system seems to be 
analogous to this latter scenario: the (R-) charge of the scalar in our case is fixed by the IIB 
construction to be 2 and it is plausible that this is a low enough value so that the transition 
remains second order all through. 

In [22], the phases of the (unbackreacted) superfluid for various values of the masses 
of the scalar field in AdS^ were investigated and it was found that for high enough mass, 
there is always a second order transition close to the normal phase. Since the probe limit 
is a large charge limit, we should expect a similar structure also in the backreacted case 
when the charge is large. That is, when the charge and the mass are both large, we should 
expect a persistent second order transition. In our IIB case, we are exploring the opposite 
limit, namely low (R-)charge and low mass (since the charge and mass are related for chiral 
primaries). Again, we find that the second order transition exists irrespective of the velocity. 
Based on these observations, it is tempting to make the suggestion that whenever the mass 
and charge are scaled together in some appropriate way, the second order transition persists 
for all velocities. Of course, to make and/or establish a precise statement along these lines 
will require a much more thorough exploration of the masses and charges of the scalars than 
we have undertaken here. Moreover, as already noticed, the persistence of the second order 
transition was also found in the AdS± case for small charges and small mass [21], while it 
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Figure 2: Free energy plots for various velocities £ = 0,0.1,0.33,0.4,0.5 (dashed lines from 
bottom to top). The RN-AdS black hole is also presented for comparison (solid line). The 
plots show that for any velocities the phase transition is second order. The apparent overlap 
of the £ = 0.5 curve with the normal phase, is an artifact of the resolution of the figure. 

was found not to exist for any value of the mass in the probe limit |22j. So it is clear that 
the appropriate statement, if it exists, will have to be dimension-dependent. 

5 Zero Temperature Limit 

One of the advantages of having a fully backreacted solution is that one can reliably go 
to the zero temperature limit. At zero velocity, the zero temperature solution is expected 
to be described by a domain wall, corresponding to the symmetry-breaking vacuum of the 
scalar potential that restores conformal symmetry in the IR. Such domain wall solution 
was constructed in |23j, and conjectured to correspond to the ground state of the type IIB 
holographic superconductor. Since we have here fully backreacted solutions at non-zero 
velocity, a natural question one would like to answer is whether and how such IR behavior 
gets modified when the superfluid flows. 

As a warm-up, and for later comparison, let us first consider the static case. A prelim- 
inary check one can perform is to see whether for £ = our condensate value tends to the 
condensate value found in [23J. This is indeed the case: for the lowest temperature point 
(T//i = 3.05 • 10~ 4 ), our condensate in the normalizations of [23] 



is ~ 0.3215, which is close enough to the zero-temperature value of ~ 0.322 found in [23 




{0)dw 



^3 



(5.1) 



(2/VV3) 3 
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R T/ / u = 8-10- 4 , ^ = R abcd R^j JT{jm = 8 ■ 1(T 4 , £ = 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 f 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 f 



Figure 3: Ricci scalar R and R a b c dR abcd as a function of the radial coordinate near the 
horizon, at zero superfluid velocity. The horizontal dashed lines mark the corresponding 
values of R and R a bcdR abcd for the UV and IR AdS geometries. 



Even without explicitly constructing the domain wall solution, one can find evidence for 
its existence by investigating the horizon values of the curvature scalars R and R a bcdR abcd - 
This strategy was adopted in [12J for superconductors in M-theory, and it was found that 
these curvature scalars on the horizon go to the AdS± values expected from a domain wall 
solution with a symmetry-breaking minimum in the IR. We can do the same computation 
here, and we do find evidence that the solution has an emergent AdS§ in the IR with the 
correct length scale. Note that the IR AdS scale, as determined by the symmetry-breaking 
vacuum [23] is L' = where we have set L — 1 in the UV. Using the fact that the Ricci 
scalar for AdS§ is — 20/ L 2 , we find that the predicted value is —22.5 in the IR. A similar 
computation using the R a bcdR abcd — 40/L 4 shows that in the zero temperature limit we 
should get the value 50.625. We plot the results for both curvature scalars in figure 3. The 
plots clearly demonstrate that at low temperatures the curvatures indeed stabilize to the 
expected domain wall values in the infrared. 

The behavior of R a b c dR abcd deserves a closer look, however. A distinctive feature of the 
present five-dimensional compared to the four-dimensional model of [12J, is that 

RabcdR abcd stabilizes to the domain wall value close to the horizon, but it starts increasing 
as the radius is further reduced. At the horizon its value is (of course) finite, but is well 
on its way to the divergence at the singularity inside the horizon! Note that in order to 

2 This sharp ascent in the curvature scalars close to the horizon is not a peculiarity of the broken phase: it 
is also there in the normal phase. For instance, R a bcdR abcd , whose expression for the normal phase Reissner- 
Nordstorm black hole we report in eq. (|D.4[) , has a similar sharp ascent at the horizon, while remaining finite 
there. On the other hand, the AdS^ case is somewhat special in that the Ricci scalar is a constant in the 
normal phase due to the tracelessness of the electromagnetic stress tensor in four dimensions. 
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make the connection with the domain wall, what we really need is the emergence of an AdS§ 
throat of the correct length scale at zero temperature, and our plots give evidence for that. 
Figure 4 reports the behavior of R a bcdR abcd zooming in near the horizon region for different 
temperatures. Happily, as the temperature is lowered the stabilized region of the plot gets 
closer to the horizon and asymptotes to the expected AdS 5 value of 50.625. 

T} jjabcd 
ttabcdft 



51.8 






-T|— 1 1 | ■ 1 1 | 1 1 1 | 1 1 1 | 1 1 I | 1 1 1 1 1 
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1 ! 


50.6 
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n 

Figure 4: Behavior of R a bcdR ahcd near the horizon for different temperatures, from left to 
right: T/fi = 3.05 • 10~ 4 (black), 1.55 • 10~ 3 (red), 3.33 • 10~ 3 (blue). The dashed vertical 
lines correspond to the corresponding horizon radii. The stabilized (i.e. domain wall) value 
RabcdR abcd = 50.625 is indicated by the dashed horizontal line. 

Let us now consider the cases with velocity, £ 7^ 0. We report in figure 5 the plot for 
the Ricci scalar vs. radius for different superfluid velocities (including the zero- velocity case, 
to ease the comparison) and in figure 6 that for R a b c dR abcd . The presence of a new scale 
means that there is a possibility that the emergent conformal symmetry in the IR is broken. 
While for low velocities our plots suggest that the same IR fixed point as the static case 
is recovered, interestingly enough, we find that for high enough velocities the conformal 
symmetry of the solution is indeed broken and the curvature scalars diverge without any 
stabilization whatsoever. This is analogous to the phenomenological models with no quantum 
critical point in the IR. The conclusion seems to be that the solutions do not stabilize to the 
conformal quantum critical point when the velocity is high enough. 

While we have not performed an exhaustive scan of velocities in this paper, it would 
be interesting to see for what precise value of the velocity this qualitative change happens, 
and study the precise nature of the phases and phase transitions, if any, there. From our 
analysis, it appears that the regime where this transition happens is between £ = 0.33 and 
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T//j = 8 ■ 1CT 4 , £ = 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 f 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 f 




R t//j = 8 ■ i(r 4 , { = 0.10 




1.0 1.2 1.4 J* 




0.2 0.4 0.6 0. 



Figure 5: Ricci scalar R as a function of the radial coordinate near the horizon for a low 
temperature. The horizontal dashed lines indicate the corresponding values of R for the UV 
and IR AdS geometries. 



Rg b c d R abcd T/n = 8 ■ IP" 4 , g = R abcd R abcd T/fi = &- 10~ 4 , j = 0.10 
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H A 4 

R„^ R abcd . r/ M = 8 ; lQ- 4 { = 0.33 R abcd R ahcd r/A* = 8 ■ HT 4 , { = 0.40 
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Figure 6: R a bcdR abcd as a function of the radial coordinate near the horizon for a low tem- 
perature. The horizontal dashed lines indicate the corresponding values of R a b c dR abcd for the 
UV and IR AdS geometries. The stabilization to the IR value, when it happens, holds till 
very close to the horizon. 
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Figure 7: Plots of the superfluid fraction vs. temperature for various values of the velocity 
f = 0.1, 0.33, 0.4, 0.5 (from right to left). 

£ = 0.4. Related to this is the observation that the condensate 



for small enough velocity. This is again indicative of a quantum phase transition: there is 
a change in the nature of the solution as we tune an order parameter at zero temperature. 
The results we find are consistent with the idea that the phase structure in the temperature- 
velocity plane is determined by the quantum critical point. It is intriguing that the relevant 
condensate seems to be measured in units of chemical potential as seen in a frame comoving 
with the superfluid flow. For a timelike vector, which for us is the superfluid velocity 4- vector, 
the time component in the rest (i.e., comoving) frame is nothing but its norm. Therefore, 
since we want to plot a scalar quantity for the dimensionless condensate, this is the natural 
choice. But unlike in the case of an ordinary fluid where the fluid velocity can be interpreted 
as arising from a boost of a static black hole, here the anisotropic part of the metric does 
not seem to have such a simple interpretation in the bulk. We intend to come back to some 
of these questions in the near future. 

Another quantity of interest^ in understanding the zero temperature limit is the superfluid 
fraction (. It corresponds to the ratio between the charge density of the superfluid flow and 
the total charge density of the system. In appendix [UJ following [2T], we elaborate on 
the interpretation of the boundary theory in terms of a two-fluid model and compute the 
expression of the superfluid fraction in terms of the fall-offs of the bulk fields, eqs. (13. 18ft - 

3 We thank Julian Sonner for raising this point. 



(0)1/3 



(5.2) 




that we plotted earlier, tends to the same value at the horizon for all values of the velocity, 
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lpT2I) . The result is 



c 




(5.3) 



This quantity is interesting because from the curves in figures 1 and 2 of [2TJ we see that 
for the AdSi case its behavior near zero temperature captures some interesting aspects of 
the nature of the phase transitions. More specifically, together with our results in this paper 
(see figure 7), we are lead to conjecture that £ — > 1 at zero temperature for all velocities 
where the rescaled condensate value at zero temperature tends to its value at zero velocity. 
From the evidence presented in [21] one could think that the zero temperature limit of the 
superfluid fraction is correlated with the existence or not of a first order phase transition 
at high enough velocity. However, in our case we have an explicit situation where we see 
a consistently second order phase transition where the limiting value of the condensate at 
zero temperature changes qualitatively as we tune the velocity. Remarkably, we find that 



corresponding value at zero velocity. Since this condensate value captures the existence or 
not of the (anisotropic) domain wall, the natural conjecture is that ( = 1, for the domain 
wall when it exists. Notice that £ — > 1 is what one would expect for the ground state 
of a superfluid flow. What we have basically demonstrated then is that three quantities 
(namely the curvature scalar (s), the rescaled condensate and the superfluid fraction) undergo 
a qualitative change at the same velocity, as we tune the velocity. We believe this is strong 
evidence for the existence/non-existence of the domain wall as we go through that velocity. 

Despite the evidence we have presented, it should be borne in mind that the preservation 
of the conformal symmetry for low velocities is not fully established. Unlike in the zero 
velocity domain wall examples discussed in the literature, we have not constructed an explicit 
solution that has emergent conformal symmetry in the IR in the cases with (low) velocity. 
However, the fact that the curvature scalars and the condensate (15. 2p stabilize to their 
respective zero velocity values (within our numerical precision), is an indication that this 
might indeed be the case. One another caveat that we emphasize here is that the perturbative 
stability of these consistent truncations in the zero temperature limit is not settled. In 
particular, when the Sasaki-Einstein manifold is a sphere, instabilities are known to exist in 
the zero temperature domain wal^lutionSEsB It is possible that for a more complicated 
choice of Sasaki-Einstein space (which is indeed what we need to have here anyway, in order 
to let the scalar chiral primary we focus on to be the operator responsible for the black 
hole phase transition (TTJ) the five-dimensional theory is stable. It is also interesting that 
the simple stringy consistent truncations do give rise to scalar potentials with symmetry 
breaking vacua, resulting in an emergent conformal symmetry in the IR at zero temperature. 

4 A related instability was recently shown to exist also in M-theory |27j for a similar consistent truncation 
for the ground state of a 2+1 dimensional superconducting system [2"51 125] . 



£ — > 1, only in those cases where the zero temperature condensate value 



takes its 
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This is precisely what one expects in the zero temperature limit of a high-T c superconductor, 
which is believed to be governed by a quantum critical point. So our expectation is that in 
the (unlikely?) event that no Sasaki-Einstein truncation can be made stable, these models 
should still capture some generic features of a holographic superffuid with emergent conformal 
symmetry in the im 
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A Asymptotic Relations 

In this appendix, we present the relations defining the dependent coefficients in the 
asymptotic expansion in the IIB case 



U = -r^o f96C C 4 /io + 48B 4 h 4 + 96C 2 h h 4 + 96B h 3 h A + C 2 h 2 ^ + B h&t + 
4oCq V 

+2AC 2 h 2 ^ + 24B hfy 1 iP 3 - Ylh\l}^\Alp + 24C h 2 L^iA xfi A tfl + ^o^ViA'o) > 

(A.l) 

fi = ( - C 2 h&t ~ Boh&t + W,hlL^lA 2 tfi + 

+Cl{hli,t - 3LVi 2 A 2 j0 ) + BohZih&l - 3LV?Xjo)) , (A.2) 

h = ~^f, h\ = kWl , Bi = £VMio, Cl = -LV?>Wb, (A.3) 

, 2(C ^? + £ ^i 3 + 3/^^ 
^3 = o^ 2 I d ^ > ( A - 4 ) 



3(C 2 + B h 2 ) 



, __3ijfA h o A i 3gA^ 

■fit 1 — „ i A„ o 



H,2 ~ 2 ' x ' 2 



(A.5) 



3 We thank Nikolay Bobev and Chris Herzog for a discussion on this point. 
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These are the general expressions when ipi ^ 0. Our primary interest will be to shoot for 
the case ipi = 0, in which case all of the coefficients above vanish identically, except 

, _ faCp - 2C C 4 hl - B 4 hl 

hi ~ 2h (C* + B hl) ■ ^ 

In particular, all the logarithmic terms vanish and we end up with a usual asymptotic 
expansion in 1/r, as expected. Note also that in asymptotically AdS solutions, Cq = as 
well. Moreover, when there is no superfluid velocity and the isotropy is not broken, B± = 
and therefore we end up getting = 0. This last result is useful in making comparisons 
with the holographic superconductor case investigated in [TTJ. 

B On-Shell Action and Counter-Terms 

In order to compute the free energy, we need the on-shell action for the type IIB system. 
As we show below it turns out that, remarkably, the on-shell action can be written purely 
as a boundary piece, and be easily evaluated. However, this boundary term is divergent: to 
cancel it we need to introduce boundary counter-terms. In what follows, we describe both 
these steps. 

For the ansatz that we work with, it can be checked directly that, despite the complica- 
tions of the equations of motion, the following relations hold 

2L 2 2L 2 
C o~ R = — T yy = — T ^- 

Here T stands for the stress tensor arising from our IIB Lagrangian, £ is defined via 

Shb = J d*xV=jj£ , (B.2) 

and R is the Ricci scalar. Notice that these relations only depend on our ansatz, i.e. they 
are true before we use the equations of motion. Going on-shell, we replace T yy and T zz by 
E yy and E zz , where E denotes the Einstein tensor E a b = R a b — \g a bR- Together with the 
relation 

E\ = -\R (B.3) 
that is valid in five dimensions, this implies that 



-g£> 







- g (^(E yy + E zz )- 2 -E a )j . (B.4) 



The right-hand-side depends only on the metric functions and can be evaluated explicitly 
for our ansatz. Direct computation reveals that it can be written as a total differential so 
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that the (on-shell) action takes the form 



,V/ // ,os = -vol i rdr(^l^V=9) , where y/=g = ^ J^- + Dir) . (D.O) 



VL 2 /i(r) 2V V ' v ^ L 3 y f{r) 

and the prime denotes the derivative with respect to r. Because of the presence of /, this 
expression is zero at the horizon and that end of the integral is safe. But it clearly gets 
contributions from the boundary, where it diverges as r 4 and we need to regulate it with 
appropriate counter-terms. 

The counter-termJ§ for the gravitational part of the action in asymptotically AdS spaces 
can be looked up in [29] . Along with these we also have to add counter-terms for the scalar 
part. The final form of these terms in our notations and conventions can be written as 



3 
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= 2 J <txy/=j\K: - -J + J Av^y- • (B.6) 

The sign convention for the extrinsic curvature is chosen so that with the outward pointing 
normal n a , 

K ab = ^(V a n b + V b n a ). (B.7) 

Note that the general gravitational counter-term discussed in [29] involves a boundary Ricci 
scalar as well: but this does not contribute for us, because our boundary becomes flat as we 
take it to infinity. The various quantities (including the scalar extrinsic curvature) can be 
computed by cutting off the spacetime at some finite r = r , then taking the limit r — > oo 
for the quantity Shb,os + Set at the end of the computation. If we define the boundary at 
r = r , then the outward normal to the surface $(£, r,x,y,z) = r — r = is n a ~ V a $, and 
after normalizing so that g ab n a n b = 1, we get 

n a = ( 0,-^L, 0,0,0 | . (B.8) 

V r VfF) J 

Since we need only the scalar extrinsic curvature, we don't need to introduce 4-D coordinates 
on the boundary and can compute it directly in the bulk coordinates as 

K ^ (8C2 ± rfB ' ± 2rC °' ± 8Bf ± rBr) (TK Ql 

K = 9 VaUh = 2L(C* + Bf)h • (R9) 

So the final form of the counter-term action is 



Set = VoL lim 



r 4 fV 2 (8C 2 + rfB' + 2rCC' + 8Bf + rBf) r 4 r— — /6 tp 2 



L 5 h^/C 2 + Bf £ 4 J \L L 



(B.10) 



6 We loosely refer to the Gibbons-Hawking term also as a counter-term, even though strictly speaking it 
is a boundary term necessary to make the variational problem well-defined. 
Note that the boundary is timelike, so it has a spacelike normal. 
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With the addition of this piece, the renormalized action Sub,os + Set no longer has the r 4 
divergence and is finite. The net result is 



S ren = vol 4 lim 



iyi/2 ( gC 2 + r j B , + 2r cC + 8Bf + rBf) 
L 5 h^/C 2 + Bj 



+ B(r) 



(B.11) 



It is interesting to note that since we are always working with solutions with ipi = 0, the 
scalar piece can in fact be omitted if one desires. 



C Superfluid Fraction 



In this section we present some details of the definition and computation of the super- 
fluid fraction ( for our solutions. We start with the renormalized action from the previous 
appendix and compute the boundary stress tensor and the boundary current by varying with 
respect to the boundary metric and the boundary components of the vector potential. 



T 



5S 



7 £7 



Ju 



6S 



(C.l) 



where now S = Sub + S ct with Sub defined by (12.11) and S ct defined by (IB. 61) . In particular, 
the relations above are not tied to our ansatz. To compute the boundary stress tensor and 
current, we need to introduce coordinates on the boundary, and we will use Greek indices 
for them. After doing the variations, using our ansatz and going on shell on the bulk, the 
resulting stress tensor and current vanish in the strict r — > 00 limit. This is consistent with 
the fact that they should be finite since we are using the renormalized action to compute 
them. The more interesting quantity is the boundary fluid stress tensor and the fluid current, 
which are defined in AdS^ via 



T, 



lim r % v , 

r— >oo 



J u = lim r 2 J u 



(C.2) 



We are using units where 167rG = 1 = L in this section. Suppressing the details and 
restricting to our ansatz these quantities can be explicitly computed in terms of the boundary 
fall-offs of eqs. IjPBjl - fEESI)) to be 



T, 



( 3/4 - B A 4C 4 
4C 4 U ~ 35 4 
5 4 + / 4 









V 















Ju 



( A t,2 \ 

A x , 2 


V / 



(C.3) 
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Now we follow the interpretation of [21] for these quantities in terms of a two-fluid model on 
the boundary, where one component is an ordinary (ideal) fluid and the other is a superfluid. 
First we can write these quantities suggestively in terms of = (—1,0,0,0) and n M = 
(0,1,0,0) as 

T^v = (e + P) Ufj,u u + Prj^ - AB^n^nv - SC^u^n^ , J^ = pu^- J s n^. (C.4) 
where 

P = U + B A , £ = 3/4-54, p = -^A t>2 , J s = --A x>2 . (C.5) 

Note that what we have done is merely to rewrite the expressions covariantly in terms of 
the vectors and n^. Another way to state the same thing is that (for example) the most 
general symmetric second rank tensor constructed from and will have to be a linear 
combination of rj^ u , u v , n u ^ and n u . 

The two fluid model can be defined by the stress tensor 

T^ u = (eo + P )u„u u + P r]^ + pp s v^v u , J /t = p n u^ + p s v^ , (C.6) 

where the subscripts n and s stand for the normal and superfluid components of the charge 
density, with the total charge density p = p s + p n . Aside from the various thermodynamical 
state variables (whose precise interpretations will not be important to us, see [21]), we have 
also introduced the superfluid velocity v ^ that satisfies the constraint ( "Josephson equation" ) 

u*X = -l. (C.7) 

The superfluid fraction is defined as 

C = — . (C.8) 
P 

Our stress tensor (IC.4p can be brought to the two-fluid form by defining t> M as 

^ = w M + 7T r V ( c - 9 ) 
o 4 

This automatically satisfies v^u^ = —1 as a consequence of = —1, and n M M M = 0. 

Rewriting our stress and current tensors (1C.4j) in these new variables we get the two-fluid 
form (ICT : 

= (e + P + 4 C 2 JB 4 ) u^u u + P % „ - (4 C 2 J B A ) v v , (CIO) 

J» = (p + J s C 4 /B 4 ) u„ - ( J s C 4 /B 4 ) Vlt . (C.ll) 

Reading off the superfluid fraction from this, we find that 

_ -{J S C 4 / B 4 ) _ A Xt2 C 4 
p A t ,2 B 4 ' 

where we have written the final result in terms of the fall-offs obtained directly from the 
solutions. This is the form we use for making the plots in figure 7. 
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D The Hairless Solution: Reissner-Nordstrom 

In understanding the phase structure, it is important to keep in mind that we are in- 
terested in comparing the free energy of the hairy black hole solution to that of Reissner- 
Nordstrom. In the five dimensional IIB case, the Reissner-Nordstrom metric [TT] can be 
given in terms of our ansatz ( 13 JJ by 

/W = l-i(l + ^)+^, A, = p(l-±), (D.1) 
h = l, B = l, C = Q, A x = 0, ip = 0. (D.2) 

In this notation, the curvature invariants studied in section [5] take the form 

R=-20f-r(10f' + rf") , (D.3) 
R abcd R abcd = 40/ 2 + 4r / (10/' + r /") + r 2 [22(/') 2 + 8r /' /" + r\f"f] . (D.4) 

All these expressions are obtained after all the necessary rescalings: we have set 167rG = 
L — th — 1- The Hawking temperature now takes the form 

r„ = i^£, (d.5) 

71 

as can be determined by the periodicity of the Euclidean section. The renormalized on-shell 
action that we determined before takes a simple form for this solution: 

S ren = 1 + 4/i 2 /9 . (D.6) 

We will compare the free energies of the hairy and hairless cases at the same T/fi to determine 
which one is the favored phase. 
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